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We postulate the existence of universal crossover functions connecting the universal parts of 
the entanglement entropy to the low temperature thermal entropy in gapless quantum many-body 
systems. These scaling functions encode the intuition that the same low energy degrees of freedom 
which control low temperature thermal physics are also responsible for the long range entanglement 
in the quantum ground state. We demonstrate the correctness of the proposed scaling form and 
determine the scaling function for certain classes of gapless systems whose low energy physics is 
described by a conformal field theory. We also use our crossover formalism to argue that local 
systems which are "natural" can violate the boundary law at most logarithmically. In particular, 
we show that several non- Fermi liquid phases of matter have entanglement entropy that is at most of 
order log (L) for a region of linear size L thereby confirming various earlier suggestions in the 
literature. We also briefly apply our crossover formalism to the study of fluctuations in conserved 
quantities and discuss some subtleties that occur in systems that spontaneously break a continuous 
symmetry. 



I. INTRODUCTION 

Recently, an exchange of ideas between quantum in- 
formation science and many-body physics has led to an 
improved understanding of the "corner" of Hilbert space 
in which ground states of local Hamiltonians reside. One 
of the most important tools for investigating the proper- 
ties of many-body ground states is entanglement entropy, 
defined as the von Neumann entropy of the reduced den- 
sity matrix of a spatial subsystem. The ubiquitous pres- 
ence of a boundary law for the entanglement entropy, 
as reviewed in Ref. has provided a rough guide to 
the entanglement properties of quantum ground states. 
This rough intuition led to a new class of quantum states 
generically called tensor network states Refs. as well 
as new insights into the classification and identification of 
many-body phases and phase transitions in Refs. l^-ITol 

In this paper we show that by considering the relation- 
ship between thermal and entanglement entropy we can 
place significant constraints on ground state entangle- 
ment structure for "natural" systems. One of our main 
motivations is to characterize the possible violations of 
the boundary law for entanglement entropy at zero tem- 
perature. There have been many constructions of anoma- 
lously entangled ground states in the quantum informa- 
tion community e.g. Ref. |Tl|, but what do these have to 
do with ordinary quantum systems relevant for labora- 
tory studies? There are also motivations from the study 
of mutual information in quantum systems at finite tem- 
perature in Refs. [T^fl3l Interesting critical phenomena 
are visible in the mutual information, and a first step 
toward understanding these numerical results is a more 
complete understanding of the temperature dependence 
of the von Neumann entropy of a single region. 

Our basic assumption that connects thermodynamics 
with entanglement is that the same low energy degrees 
of freedom are responsible for both long range entangle- 



ment and low temperature thermal physics. To give a 
concrete example, in one dimensional relativistic criti- 
cal systems, while the high energy physics contributes an 
area law term to the entanglement entropy, only the low 
energy modes contribute to the log L entanglement and 
low energy thermal properties as discussed in RefsJ^— . 
Indeed, there is a universal crossover function that in- 
terpolates between the zero temperature entanglement 
entropy and the finite temperature thermal entropy of 
a given subregion (consisting of a single interval). We 
study and generalize this crossover phenomenon in vari- 
ety of critical systems in different dimensions. 

In more detail, we will make the following assump- 
tions throughout this paper. We always study gapless 
systems since it is obvious (though not rigorous) that 
generic gapped phases obey a boundary law. Our pri- 
mary assumption is that there exists a universal crossover 
function that relates thermal and entanglement entropy. 
This crossover function is only defined up to boundary 
law terms coming from high energy physics. We also as- 
sume that the system does not possess extensive ground 
state degeneracy and that the Hamiltonian is not fine 
tuned (beyond the tuning necessary to reach criticality). 
We will mostly consider translation invariant states, but 
we do discuss disordered states in Sec. VI. In short, we 
want to consider sensible gapless ground states of local 
Hamiltonians, but in an effort to be precise, we give 
the above assumptions as sufficient criteria for "sensi- 
bleness". Finally, let us note that the renormalization 
group perspective on entanglement structure permeates 
our entire discussion. 

We study the von Neumann entropy S(L,T) = 
— Tr (pl log Pl) of a real space region of linear size L 
in d spatial dimensions as a function of temperature, 
L, and region geometry. Recall that at zero temper- 
ature most gapless quantum systems in d > 1 dimen- 
sions satisfy a boundary law for the entanglement entropy 
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Sl ~ L d 1 with the coefficient of this term non-universal 
(see Ref. Q). However, there are gapless systems that 
violate the boundary law for entanglement entropy in- 
cluding free fermions with a Fermi surface^—, Landau 
Fermi liquids 20 , and Weyl fermions in a magnetic field 
at weak and strong coupling 2 ^. These examples have an 
entanglement entropy ~ L d ~ l log (L). 

It is of enormous interest to generalize this result to 
understand the entanglement structure of non-Fermi liq- 
uid ground states of matter. Many such states share with 
the Fermi liquid the crucial feature that there are gapless 
excitations that reside at a surface in momentum space. 
However unlike in a Fermi liquid there is no description of 
these excitations in terms of a Landau quasiparticle pic- 
ture. Such states were suggested to also violate the area 
law for the entanglement entropy based on a heuristic 
argument that views that gapless momentum space sur- 
face as a collection of effective one dimensional systems™. 
If the area law is indeed violated can the violation be 
stronger than in a Fermi liquid? 

An example of a non-Fermi liquid state with a gapless 
surface in momentum space was studied numerically in 
ReL— . The second Renyi entropy of a wavefunction (ob- 
tained by Gutzwiller projecting a free Fermi sea) for a 
gapless quantum spin liquid phase of an insulating spin 
system in two dimensions was calculated using Monte 
Carlo methods. The second Renyi entropy was shown to 
obey a behavior consistent with L log (L). Given the cur- 
rent limitations on system size it is hard to distinguish 
this from a power law violation of the area law. It is 
therefore important to have a general understanding of 
how seriously the area law can be violated in such a spin 
liquid state. 

The quantum spin liquid phase discussed above is ex- 
pected to be described by a low energy effective theory 
with a Fermi surface of emergent fermionic spin- 1/2 par- 
ticles (spinons) coupled to an emergent U(l) gauge field. 
Similar effective field theories describe Bose metals, some 
quantum critical points in metals, and other exotic gap- 
less systems. In all these cases the violation of the bound- 
ary law is suggested by heuristic arguments. Based on 
the analogy with Fermi liquids we might guess that the 
violation is logarithmic^. It is clearly important to have 
a firm argument for the correctness of this guess. Provid- 
ing such an argument is one of the purposes of this paper. 
What about other non-fermi liquid states where the ef- 
fective theory is not yet understood? We will address 
a class of such states that have a critical Fermi surface 
with appropriate scaling properties^ to discuss the scal- 
ing constraints on their entanglement structure. In all 
these cases we argue that the L d_1 log (L) is the fastest 
possible parametric scaling with L in d dimensions. 

Besides the von Neumann entropy, we also investi- 
gate the scaling behavior of fluctuations in conserved 
quantities as in Refs. [THHi Here the structure is 
slightly richer, but the basic conclusions are very simi- 
lar. In phases with unbroken symmetry, the fluctuations 
in the conserved quantity generating the symmetry scale 



no faster than L d_1 log (L) at zero temperature, again 
under the assumption that the same low energy modes 
responsible for thermal fluctuations also give rise to these 
zero temperature fluctuations. 

This paper is organized as follows. We begin with a 
discussion of the crossover behavior in the simplest con- 
formally invariant case in d dimensions. Next we discuss 
the case of codimension one critical manifolds relevant 
for Fermi liquids, and then we discuss the general struc- 
ture including higher codimension critical manifolds. Fi- 
nally, we turn to a discussion of fluctuations in conserved 
quantities. We conclude with a discussion of possible vi- 
olations of our scaling formalism. 

II. SCALING FORMALISM: INTRODUCTION 
A. Conformal symmetry 

Consider a local quantum system with Hamiltonian H 
at finite temperature T = so that the entire system 
is in the mixed state p(T) cx exp (—(3H). As (3 — > we 
recover the ground state up to corrections exponential in 
the gap to the first excited state. We will be exclusively 
interested in systems where H is either in a gapless phase 
or at a critical point. Thus we will always have some 
notion of scaling symmetry although we will often not 
have the full power of the conformal group. 

Consider now a region R of linear size L inside a larger 
many-body system. The complement of region R is de- 
noted R. The reduced density matrix of R is 

p R (L,T)=Tr R (p(T)) (1) 

and the von Neumann entropy of this reduced density 
matrix is 

S R (L,T) = -Tr R (p R logp R ). (2) 

We will also be interested in generalizations of the von 
Neumann entropy called Renyi entropies labeled by a pa- 
rameter n: 

Sn = -^—Tr R (p R ). (3) 
1 — n 

The limit n — > 1 of S n is simply S R , the von Neumann 
entropy of p R . 

Let us initially consider the special case of a conformal 
field theory in d spatial dimensions. Two simple limits ex- 
ist. As LT — > (a non-universal velocity is suppressed) 
the von Neumann entropy recovers the usual entangle- 
ment entropy of the ground state. We know from ear- 
lier studies that the entanglement entropy may contain 
a mixture of universal and cutoff dependent terms, see 
Refs. I26I4281 for representative calculations. For example, 
the boundary law term, going as L 1 is non- universal, 
but there are universal logarithmic terms in d = 1, 3, 5, ... 
dimensions. In d — 2, 4, ... dimensions the universal log- 
arithmic term is replaced by a constant term. On the 
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other hand, as LT — >• oo we recover the usual thermal 
entropy going as (LT) d . 

Using our basic assumption we write the entropy of 
region R as 

S R {L,T)=T*f R {LT) + ... (4) 

where ... stands for the aforementioned addition or sub- 
traction of non-universal terms involving the momentum 
cutoff A. Let us now determine the properties of Jr and 
the exponent (f>. For the moment we suppress the re- 
gion dependence writing fa as /. First, as LT ->oowe 
must recover the extensive thermal entropy and hence 
f(x — > oo) ~ x d . This further implies that = to ob- 
tain the correct temperature dependence of the entropy. 
In the opposite limit as LT — > the only possibility for 
a non-zero and finite contribution is f(x) — > constant or 
fix) — > \ogx. The possibility of the logarithm is allowed 
since the T appearing in the logarithm can be replaced 
by A at the expense of a non-universal term. 

We conclude that our scaling assumption is consistent 
with either a universal constant term or a universal log- 
arithm in the entanglement entropy of a conformal field 
theory at zero temperature. Indeed, both these possibil- 
ities are realized, the logarithmic term obtains for d odd 
and the constant for d even. This is also an appropriate 
time to mention the possibility of shape dependence, for 
example, the fact that sharp corners produce logarith- 
mic corrections is completely consistent with our scaling 
framework (see Refs. l27l[29h . It is important to note that 
the constant term in odd dimensions is only meaningful 
in the absence of corners. Other types of universal terms 
like (logL)P ^ 1) or L d - 1+S (S > 0) are not allowed 
unless they violate our assumptions and are unrelated to 
the thermal physics. Of course, this conclusion is very 
natural from the renormalization group point of view. 

We briefly elaborate on this point and discuss the 
structure of high energy contributions in more detail. 
Locality demands that all high energy contributions be 
proportional to integrals of local geometric data over the 
boundary. Consider an entangling surface OR in d di- 
mensions. We may use coordinates u a (a = 1, ...,d — 1) 
in terms of which the surface is x l (u a ) and the induced 
metric is h a b = d a x l dbX J 8ij (we only consider flat space 
here, the generalization is straightforward). In addition 
to the intrinsic geometry of the surface we also have ex- 
trinsic geometry related to the embedding of the surface 
into flat space. For example, a cylinder has extrinsic cur- 
vature but no intrinsic curvature where as a sphere has 
both. The extrinsic geometry is controlled by the extrin- 
sic curvature which is given in terms of the normal vector 
n l and the projector onto the surface = <5* — rfrij as 
Kij = P^dkrij. Now an important constraint for global 
pure states is the requirement that S{R) = S{R), and 
since we consider here only high energy contributions 
that are independent of the low energy physics, we may 
still demand this symmetry at finite temperature for the 
terms of interest. Since the only difference between the 
boundary of R and R is the direction of the normal n we 



conclude that only even powers of n can appear. Thus 
only even powers of the extrinsic curvature and hence 
only even powers of derivatives can appear (the same 
is true for intrinsic terms). Roughly speaking, we must 
form fully contracted invariants involving the normal vec- 
tor and the gradient, but requiring the normal vector to 
appear with only even powers forces the same for gradi- 
ents due to rotation invariance. This explains the general 
even/odd structure of universal terms via a simple scaling 
argument e.g. Ref. H(J 

Let r denote the length scale of interesting along the 
RG flow. The infinitesimal contribution to the entangle- 
ment entropy from degrees of freedom at scale r is of the 
form described above: 

dS L d - 1 ( r 2 \ 

where we have used the appropriate logarithmic mea- 
sure for r. The presence of only even corrections comes 
from our argument above. Performing this integral from 
the UV r = e to the IR r = L gives the desired struc- 
ture. More generally, one should cut this integral off at 
min (L, 1/T) where, roughly speaking, the von Neumann 
entropy becomes thermal in nature. 

Returning to our main development, the simplest ex- 
ample of such a crossover function occurs in d = 1 confor- 
mal field theories where the single interval case is dictated 
by conformal invariance. The result is 

(6) 

and we see immediately that this form is consistent with 
our general scaling hypothesis. The high energy cutoff 
A can be shifted by a boundary law respecting term, 
but the thermal physics and long range entanglement is 
independent of the precise choice of A (again up to a 
non- universal dimensionful conversion factor) . 

B. Conformal field theories in d > 1: Holographic 
calculation 

We have already mentioned one concrete example of 
such a crossover function in 1 + 1 dimensional confor- 
mal field theory. In general the computation of such a 
crossover function is a highly non-trivial task for interact- 
ing conformal field theories in spatial dimensions d > 1. 
However, one set of examples where a computation is 
possible is provided by holography. For an introduction 
to this set of ideas, see Ref. [31]. For our purposes it suf- 
fices to mention three facts. First, certain strongly inter- 
acting conformal theories living in flat space are dual to 
theories of gravity in a curved higher dimensional space 
known as Anti-de Sitter (AdS) space . The UV of the 
conformal field theory lives at the boundary of the grav- 
itational spacetime. Second, there is a simple prescrip- 
tion to compute the entropy in such theories, at least 
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in a special limit. Briefly, we must compute the "area" 
(in Planck units) of a minimal "surface" in the extended 
gravitational geometry that terminates in the UV on the 
boundary of the region of interest. Third, finite tem- 
perature effects are dual to placing a black hole in the 
gravitational spacetime. 

Putting all these facts together permits a geometric 
calculation of the entropy of the field theory that is pre- 
cisely of the form we assumed. We will now give the 
details of this calculation. Consider first the case of 
CFT 1+1 . x and t are the CFT directions while r is the 
emergent scale coordinate (r — > is the UV boundary). 
The gravitational geometry in this case is the AdS2+i 
black hole with metric 



T 2 



-fdt 2 



dr 2 + dx' 



(7) 



where f(r) = 1 — r 2 /r 2 ,. We study minimal curves r(x) 
terminating in an interval of length £ at finite tempera- 
ture. The minimal length is given by 



ft/ 2 dx , 

/ — VT+W/dxYJ^. 



(8) 



Minimizing this length with respect to r gives an equation 
of motion which may immediately be integrated to yield 
a conserved quantity 



1 



f 



ry/l + (dr/dx) 2 f- 



(9) 



where r m is the maximum depth achieved by the curve. 
Solving this for dr/dx gives 



dr/dx = 




1 / = 



We may now rewrite the length as 

dr r r„ 



J = 2 



j: 



V(r 2 m -r 2 )fr 2 
while the parameter r m is determined from 

1=2 1 !hr 



V(r 2 m ~r 2 )f 



(10) 



(11) 



(12) 



e is a UV cutoff. 

Changing variables to w = r 2 jr 2 m allows us to rewrite 
the integral for £ as 

r 1 i 

£ = r dw (13) 
Jo v(l - w)(a 2 - w) 

where we have safely put e = and a = r /r m > 1. 
Performing the integral we obtain 

/^logf^Wogf^Y (14) 



{a -If 



a-l 



We may now solve for a in terms of £ 

1 + e - £ / ro 



a = 



1 _ e -t/r ■ 

Returning now to the length J we find 

J=— [ dw — = — . 

r m Je \/(l - w)(a 2 - w) w 



Doing this integral gives 



(15) 



(16) 



(17) 



To compute the entropy we append the factor and 



use the relation c = to obtain 



5 = i l0g (v sinh (^ 



(18) 



We may now determine r$ in terms of the temperature. 
Zooming in near the horizon and writing r = ro — p we 
have 



ds 2 ^%(J± dt 2 + ^ dp 2 
rn V r 2p 



'o 



(19) 



Changing variables to u = 2^fp the near horizon metric 
is brought into the form 



ds 2 %dt 2 + du 2 . 



(20) 



Demanding periodicity in imaginary time gives (3/r — 
2tt or 2r = l/(nT). Plugging this into our entropy for- 
mula reproduces the usual crossover function 



S(£,T) = ^log(-^- S mh(7rT£) 



(21) 



Let us now consider regions in higher dimensional con- 
formal field theories. For a d + 1 dimensional CFT the 
relevant metric is the AdS^+2 black hole 



r2 

ds 2 = -£ ( -fdt 2 + - £ dr 2 + dx 2 



1 

r 



(22) 



where f(r) = 1 - r d+1 /r% +1 . 

We consider strip-like regions in the field theory of 
cross-section A and width I (A has units of length'' -1 ). 
The minimal surface area is 



ft/ 2 dr , 

u = A\ —^1 + (dr/dx) 2 f-K (23) 

J -1/2 r 



As before, we may immediately integrate the equation of 
motion to yield the conserved quantity 



1 



VI + (dr/dx) 2 /- 1 r d ' 



(24) 
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r m has the same meaning as before. Solving for dr/dx 
we find 



dr/dx 



y/(r%-r**)f 



Putting these facts together gives two integrals 
' drr d 1 r£ 



a = 2A 



V(r™-rM)fr d r d 



and 



dr r d 



y/(r%-r 2d )f 



(26) 



(27) 



We may set e = in the integral for I which gives a 
cutoff independent relation r m — r m {ro,t). Of course, 
when To — > oo(/ — > 1) we have r m oc I. 



For the area integral we write 



(25) — = 

2A 



dr 




(28) 

We have essentially subtracted off the UV sensitive 
boundary law term so that the integral in parenthesis 
converges as e — > 0. The e dependence is now trivial to 
extract and wc find 



IT 

A 



1 e 



d-l 



(a/A) & J£,r ) 



(29) 



But this is of the required cross-over form: a universal 
cross-over function plus a boundary law respecting piece 
sensitive to UV physics. 



We may put this function into the precise form considered above by scaling the variables appropriately. In (j2"7| set 
r m = tyi^T) and w — r/r m to obtain 



1 = g(£T) / dw 



o Vl - w 2d ^l - k{£Tg{£T)) d +' s -w 



d+l 



(30) 



with k some constant. This is an implicit equation for the scaling function g(x) that can easily be shown to have the 
properties claimed above. In particular, it shows relativistic length-energy scaling. Plugging this scaling form into 
(l28l) shows that (a/A)n n has the form 



1 



(£Tg(£T)) 



d-l 



dw 



1 



■j d Vl - w 2d y/l ~ k(£Tg{£T)) d + 1 w d + 1 w d 



d-l 



T d -'f{£T). (31) 



One comment is necessary, the overall factor T d-1 differs from the 4> obtained above simply because we are here 
working in a limit where A is bigger than all other scales. Repeating our general analysis above predicts <fi = d — 1 
since we must have f(x — > oo) — > x. Similarly, we have f(x — > 0) — > x - ^ -1 ) to compensate the vanishing powers of 

I 



C. Non-relativistic scale invariance 

In this section we discuss critical theories with dynami- 
cal exponent z ^ 1. The dynamical exponent controls the 
relative scaling of space and time leading to the invari- 
ant form LT 1 / 2 where again we suppress a non-universal 
dimensionful parameter. The thermal entropy of such a 
theory scales as L d T d / z as follows simply from the re- 
quirement of dimensionlessness and extensivity. Let us 
again introduce a universal scaling function following the 
assumptions above. We write the entropy as 

S(L, T) - T^/(iT 1/z ) + ... 

and use the limit LT X I Z — > oo to establish that f(x) — > 
x d and (f> — 0. The rest of the analysis for the non- 
relativistic case is unchanged and again we are permitted 
at most universal constant or logarithmic terms in the 
entanglement entropy. 



One can also perform a holographic computation in 
this setting using so called Lifshitz geometries. In fact, 
the spatial part of the metric is unchanged at zero tem- 
perature, hence the structure of the entanglement en- 
tropy is identical. For example, even space dimensions 
have subleading constants while odd space dimensions 
have subleading logs. These zero temperature solutions 
may also be generalized into finite temperature black 
holes solutions, at least for certain values of z. Ref. HU 
contains a nice example of such a Lifshitz black hole with 
d = z = 2 and metric 

,2 1 dr 2 dx\ . . 

ds =-f^I + J— + 7^ (32) 

where / = 1 — r 2 jr^. As claimed, the only difference be- 
tween this metric and the relativistic examples above is in 
the r dependence of the dt 2 term and the different power 
of r appearing in /. The same manipulations establish 
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a nearly identical crossover structure to the relativistic 
case except that the argument of all scaling functions is 
iT 1 / 2 instead of IT, As always, a dimensionful constant 
has been suppressed. 

III. SCALING FORMALISM: CO-DIMENSION 1 

Now we turn to the case where the low energy degrees 
of freedom reside on a codimension 1 subspace in mo- 
mentum space. By contrast, the scale invariant theories 
in the previous section had low energy degrees of free- 
dom only a single point in momentum space (or finite 
set of isolated points) . Examples of systems with a codi- 
mension 1 gapless surface includes Fermi liquids with a 
d — 1 dimensional Fermi surface in d dimensions, Bose 
metals, spinon Fermi surfaces, and much more. Later we 
will consider the case of a general codimension gapless 
manifold. 



A. Review of Fermi liquids 

The low energy physics of a Fermi liquid is, for many 
purposes, effectively one dimensional (an exception to 
this rule is provided by zero sound which requires the full 
Fermi surface to participate). Thus Fermi liquids violate 
the boundary law for entanglement entropy because one 
dimensional gapless systems violate the boundary law. 
The anomalous term has been found to be universal in 
the sense that it depends only on the geometry of the 
interacting Fermi surface and not on any Landau param- 
eters. Remarkably, this term also controls the finite tem- 
perature entropy to leading order in T/Ep. The universal 
part of the entanglement entropy and the low tempera- 
ture thermal entropy are connected by a universal scaling 
function which can be calculated using one dimensional 



conformal field theory. 

We work in d = 2 for concreteness. Consider a real 
space region A of linear size L in a Fermi liquid with 
spherical Fermi sea T. The entanglement entropy for 
this region scales as Sl ~ kpL log L + L/e + ... with the 
boundary law violating term universal and the subleading 
term non-universal. The precise value of the boundary 
law violating term is expressed in terms of the geometry 
of the real space boundary OA and the Fermi surface dT 
as 

Sl = ^—^[ [ dA x dA k \n x -n k \\og(L) (33) 
27r 12 JdA Jar 

where n x and n k are unit normals to d A and dT. The 
intuition behind this formula is simply that the Fermi 
surface in a box of size L is equivalent to roughly IcfL 
gapless modes that each contribute logL to the entan- 
glement entropy. This formula is known as the Widom 
formula because of its relation to a conjecture of Widom 
in signal processing. The Widom formula has not yet 
been rigorously proven, but it has been checked numer- 
ically and can be obtained simply from the one dimen- 
sional point of view. To generalize to finite temperature 
we must replace the zero temperature one dimensional 
entanglement entropy by the general result at finite tem- 
perature given by 

Sl+l(t ,^£±f log (^ sinl (-)). (34) 

Fermi liquids are described by many nearly free chiral 
fermions with c = 1 and c = 0. The marginal forward 
scattering interactions do not change the number of low 
energy modes, and hence the mode counting picture still 
works quantitatively. However, we will only require a 
much more crude scaling assumption for our results. 



Following this one dimensional result a higher dimensional Fermi liquid also possesses a universal crossover between 
the low temperature thermal entropy and the universal part of the entanglement entropy. This scaling function 
depends only on the geometry of the real space region A and on the shape of the Fermi surface dT. For a spherical 
Fermi surface and spherical real space region of radius L this universal crossover function is given by 

S(L,T) = ±±- [ [ dA x dA k \n x - n k | log ( sinh ( ^fK ' "*h ) . (35) 



27T 12 JqaJqt V V PvF(n k 



B. Non- Fermi liquids 

We now consider the entanglement structure of non- 
Fermi liquid states. We will restrict attention to the class 
of such states that have a codimension 1 gapless surface 
in momentum space, a critical Fermi surface, but where 
there is no Landau quasiparticle description of the exci- 
tations. A general scaling formalism has been developed 
for these states in Ref. [23, Of primary importance to 



us is the scaling of the thermal entropy. Our considera- 
tions will apply to any non-Fermi liquid falling into the 
general scaling formalism of Ref. [23| irrespective of the 
detailed low energy theory. However, to be concrete let 
us consider a Fermi surface coupled to a gauge field in 
d=2. 

Recently there has been a controlled calculation of the 
properties of this system in terms of the gauge field dy- 
namical critical exponent zj, and the number TV of fermion 
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flavors in Ref. following important earlier work in 
Refs. HHl The expansion parameters are e = z b - 2 
and 1/N with a controlled limit possible as N — > oo with 
eN fixed. This system was found to possess a critical 
Fermi surface. Following the intuition for Fermi liquids, 
this system will violate the boundary law for entangle- 
ment entropy because of the presence of many gapless 
one dimensional degrees of freedom. However, this sit- 
uation is not a trivial generalization of the Fermi liquid 
case because the system lacks a quasiparticle description. 

Thermodynamic quantities can be understood roughly 
in terms of many one dimensional gapless degrees of free- 
dom on the Fermi surface with a dynamical critical expo- 
nent Zf 7^ 1. The thermodynamic entropy is predicted to 
be S ~ kpT x l zs (fcf just measures the size of the Fermi 
surface). Additionally, the low energy theory is such that 
only antipodal patches of the critical Fermi surface cou- 
ple strongly to each other. With our current knowledge, 
we cannot formulate the patch theory as a truly one 
dimensional theory, nevertheless thermodynamic quan- 
tities are correctly captured. Furthermore, although the 
Fermi surface curvature must be kept in all existing for- 
mulations, this curvature enjoys a non-renormalization 
property which makes it into a kind of gauge variable: 
dispersing perpendicular to the patch normal is roughly 
like changing patches. However, we reiterate that our re- 
sults depend only on the thermodynamic entropy being 
given by S ~ kpT 1 '^. 

As usual, we write the von Neumann entropy as 

S(L, T) = kpT^fiLT 1 ^) (36) 

where z = zj is the fermion dynamical critical expo- 
nent and 4> is an exponent to be determined, kp just 
measures the size of the Fermi surface. It does not en- 
ter into the scaling argument in a non-trivial way. Now 
from thermodynamics we know that for x = LT 1 ^ — > oo 
we must have S ~ T X I Z L 2 . The L dependence requires 
that f[x) ~ x 2 as x — > oo, and the T dependence forces 
us to choose <f> = — l/z. As x — LT X / Z -> we must 
have f(x) ~ x in order to cancel the diverging T depen- 
dence. More generally, we must have f(x) — xf(x) with 
f{x) — a + blogx as before. In particular, powers of 
log are not allowed because these would produce T — >• 
divergent terms that are supposed to be finite and univer- 
sal. This demonstrates that these non- Fermi liquid states 
may violate the boundary law at most logarithmically. 

Let us also make some more detailed speculations. The 
zt = 2 critical Fermi surface actually corresponds to a 
marginal Fermi liquid, so for this theory with N flavors 
we suspect that the boundary law violating term has the 
usual Fermi liquid form 

Sl = Nx^±-[ [ dA x dA k \n x -n k \log(L). (37) 
Ztt 12 J dA J dT 

At finite e = zj - 2 we expect modifications of the prefac- 
tor due to e dependent corrections. However, it is likely 
that the geometric dependence of the integral remains 



unchanged. Indeed, the different patches in the critical 
Fermi surface decouple much more strongly than they do 
in a Fermi liquid. We also note this geometrical form has 
recently been verified in a holographic setup with log vio- 
lations of the boundary law. Depending on how precisely 
the critical Fermi surface is effectively one dimensional, 
a more structured crossover function of the form 

Sl = ^ Tl/Z °,i N ' e) [ [ dAxdA^-n^f^iLT 1 '*) 

27r 12 JdA JdT 

(38) 

may be expected. The function / £j jv would play the role 
of log (sinh (ttx)) in the Fermi liquid case. This detailed 
geometric form may be too strong a requirement in gen- 
eral, but the general scaling form in the previous para- 
graph is certainly reasonable. 

Similar scaling arguments can be made for the Renyi 
entropy S n = log (Tr(p™)). We know that the Renyi 
entropy at finite temperature has a specific relationship 
to the thermal entropy because of the simple scaling with 
T of the finite temperature free energy. The complete 
result is 



S n (T) = -S(T) (39) 

n-l l + ± 

where S(T) is the thermal entropy. This n dependence 
of the Renyi entropy actually holds for all T and L in 
the z = 1 case of a Fermi liquid. It would interesting 
to determine if this is also true for the z ^ 1 theory. 
Because the Renyi entropy is potentially much easier to 
calculate numerically and analytically we believe it is a 
useful target for future work. 



IV. SCALING FORMALISM: GENERAL 
CO-DIMENSION 

As a simple example of what we have in mind, con- 
sider a free fermion system tuned so that the "Fermi 
manifold" is of codimension q in the d dimensional mo- 
mentum space, q = 1 is the generic case, a Fermi line in 
d = 2 and a Fermi surface in d = 3. q = 2 in d = 2 and 
q = 3 in d = 3 correspond to Dirac points. The interest- 
ing case of q = 2 in d = 3 is the problem of Fermi lines 
where the zero energy locus is one dimensional in the 
three dimensional momentum space. The codimension is 
a useful parameter because it tells us the effective space 
dimension of the local excitations in momentum space. 
For example, the Fermi surface case always q = 1 indi- 
cating that the excitations are effectively moving in one 
dimension (radially). Similarly, the case q = 2 in d = 3 
corresponds to modes that move in two effective dimen- 
sions since there is no dispersion along the Fermi line. 
Just as we could calculate the entropy of Fermi surfaces 
by integrating over the contributions of one dimensional 
degrees of freedom, we can obtain the entropy of these 
higher codimension systems by integrating over the con- 
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tributions of dimension q degrees of freedom. We now 
make this explicit with a scaling argument. 

Suppose that there exists a universal scaling function 
connecting the entanglement entropy to the thermal en- 
tropy for these codimension q free Fermi systems. If we 
assume a generic bandstructure then the dispersion may 
be linearized near the Fermi manifold yielding a dynam- 
ical exponent z = 1. Thus we write the von Neumann 
entropy of region R as 

S(L, T) = k d F q T^j(LT) (40) 

where the factor of k d F ~ q accounts for the size of the Fermi 
manifold. The thermal entropy of such a system scales 
as kp~ q L d T q and hence as usual we require f(x) — > x d as 
x — > oo. This also fixes <p = q — d < 0. Requiring regular- 
ity in the limit as T — > 0, we see that f(x) — > x d ~ q f{x) 
as x — > where f(x) = a + b log x. Thus we discover that 
such systems may have a universal term proportional to 
(kFL) d ~ q with either a constant or logarithmic prefac- 
tor. We emphasize that this is precisely what one ex- 
pects from integrating a q dimensional contribution over 
the Fermi manifold. In particular, we expect a constant 
prefactor for q even and a logarithmic prefactor for q odd 
because the q dimensional system has z = 1 and hence re- 
sembles a relativistic scale invariant theory of the type we 
considered in Sec. II. These statements may be checked 
in the free fermion case because the entanglement en- 
tropy can be computed exactly, however, we defer a full 
discussion of this case to a future publication. 

We conclude this section by noting that our conclusion 
is unmodified even if we have an interacting theory with 
a codimension q gapless manifold and with general z ^ 1. 
The scaling function has the form 

S = kjr q T*f(LT 1 /*) (41) 

with (f> = (q — d)/z which is obtained by matching to 
the thermal entropy k F ~ q L d T q / z . Although we do not 
rigorously prove that the thermal entropy is always of 
this form, such a form does follow from a very general 
scaling analysis in momentum space and we know of no 
exceptions. We still predict a universal term, constant or 
logarithmic, proportional to (fcirL) d_9 at zero tempera- 
ture. For an example of such a transition, see Ref. [36| 
which considered a quantum critical point between a line 
nodal metal and a paired superconductor. 

V. FLUCTUATIONS OF CONSERVED 
QUANTITIES 

In this section we give a brief description of our scal- 
ing formalism as applied to an interesting observable: the 



fluctuations of a conserved charge. These considerations 
are motivated by the direct experimental accessibility of 
such fluctuations as well as by a desire to illustrate the 
general nature of our arguments. As the primary ex- 
ample, consider a conserved number operator N that 
may be written as a sum of local densities N = Y] n r . 
This operator commutes with the Hamiltonian and we 
may label energy eigenstates with different values of N. 
Hence N itself need have no fluctuations in the ground 
state. However, we can consider the restricted opera- 
tor Nji = J2 r eB Ur wri ich need not commute with the 
Hamiltonian and may have fluctuations. An interesting 
measure of the correlations between R and its environ- 
ment is thus the quadratic fluctuations in Nr 



AN R = ((N R -(N R )) 2 ). (42) 



We expect in a gapped phase of matter that these fluc- 
tuations satisfy a boundary law AN R ~ L d ~ l . Gap- 
less phases in higher dimensions also appear to satisfy 
a boundary law so long as the symmetry generated by 
N is unbroken. As usual, one dimension is an excep- 
tion where it is known that AiVj| ~ logL although 
the prefactor is not as universal as for the entangle- 
ment entropy i.e. it depends on the Luttinger param- 
eter. It has also been shown that Fermi liquids violate 
the boundary law for fluctuations in higher dimensions 
with AN% ~ ^{kpLY- 1 log (L). 

Because the discussion is so similar to the case of the 
entropy, we will not give any of the details here. How- 
ever, a few points are worth mentioning. First, what we 
study is now the crossover between the zero temperature 
number fluctuations and the finite temperature number 
fluctuations as controlled by the thermodynamic com- 
pressibility. Indeed, precisely such a crossover argument 
was used in Ref. to argue that the number fluctuations 
in a Fermi liquid are modified by Fermi liquid parameters 
(unlike the entanglement entropy). As before, we are al- 
lowed to subtract off any boundary law contribution since 
they can be generated by high energy degrees of freedom. 
For relativistic scale invariant systems the finite temper- 
ature fluctuations go like (LT) d and we again find that 
only constant or logarithmic universal terms are allowed 
at T = 0. More generally, we conclude that L d ~ Y logi is 
the most severe violation of the boundary law that is pos- 
sible without breaking a symmetry e.g. in critical Fermi 
surface systems that remain compressible (see Ref. l23h . 
If the symmetry is broken then we expect AN R ~ L d as 
follows from a mean-field wave function for a superfiuid 
of the form JT e ab ^ |vac). 



For conformal field theories these statements may be straightforwardly checked using the fact that the conserved 
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current J M has dimension A = d. To compute the charge fluctuations we must compute 

AN%= [ d d x [ d d y(J°(x)J a (y)} = [ d d x [ d d y- 3_ + contact terms. (43) 

Jr Jr Jr Jr \ x ~ y\ 

The contact terms are necessary to cancel a naive short distance volume scaling in the fluctuations. Such a volume 
scaling from UV modes is unphysical because UV fluctuations are local and hence can only change the charge in 
region R by appearing near the boundary. Let us take R to be a sphere of radius L. We can do the integral over y 
by changing variables to w = y — x and integrating over w in polar coordinates. The limits of integration are from 

r = \w\ = to r = r*(#) = — \x\ cos 9 + \J L 2 — \x\ 2 sin 2 9 where 9 is the angle between w and x. Let us specialize to 
d = 2 for simplicity. The fluctuations now go like 



L**I"^ (14) 



where we have discarded the aforementioned UV divergent volume piece, r* remains non-zero for all 9 unless |x| — > L, 
and it can be shown that the leading singularity is 



/ 



i.e. the boundary law. If we further expand the 9 integral in powers of L — \x\ we also find a sub-leading term giving 
rise to a logarithmic correction in agreement with our scalin g result. 



We wish to mention one final subtle point that arises in 
phases with broken symmetry and which is not properly 
captured in our thermodynamic treatment. Thus con- 
sider a superfluid phase where the particle number sym- 
metry is broken. Besides the usual sounds modes that 
possess an energy scaling as l/£ (I is system size), there 
is also a zero mode with energy levels scaling like l/£ d . 
This zero mode plays an important role in finite size sys- 
tems, like quantum Monte Carlo simulations, where it is 
responsible for insuring that although the symmetry is 
broken, the many-body state has a definite particle num- 
ber. In other words, it is related to the fact that the 
many-body ground state in a finite size system is a prop- 
erly a cat state unless the symmetry is broken explicitly. 
This zero mode is not easily visible in thermodynamics 
but it does affect the entanglement entropy and number 
fluctuations in an important way. 

Let us ignore the sound modes and ask for a state that 
capture the dynamics of the zero mode. Such a cat state 
has the form 



\M) 



I 



d9e 



-iMO 



\0)r 



(46) 



where \9) r is a state of definite phase on site r. We wish 
to trace out part of the system and compute the entropy 
of the remainder, but this problem has to be regulated 
because ambiguities arc encountered in this procedure. 
Consider a simpler system consisting p states per site 
with a Z p symmetry relating them and where the many 
body state is of the form J2x=i ®r \ x )r- If we now trace 
out part of the system we find a reduced density matrix 



for region R of the form 



PR = toRP = ^®\x){x[ 



(47) 



x reR 



Perfect correlation with the environment has rendered 
the reduced density matrix completely diagonal. The en- 
tropy is now trivially S = log p. To connect this model to 
the superfluid, we need only estimate the effective value 
of p. We do this by counting the effective number of or- 
thogonal states in R. Now the many-body coherent state 
of the form 



ae 



i^) = (g) e - |Q|/2 E : ^i^ 



r£R 



(48) 



has an overlap with a neighboring state of the form 



(9\9') = cxp[\R\\a\ 2 (e-* (e - e "> - 



(49) 



Expanding in small 9 — 6', the first real term is 
exp — \R\\a\ 2 (9 — 9') 2 /2 and hence states greater than 
A9 - 1/vTRIH 2 are effectively orthogonal. Hence we 
may take p ~ 2n/A9 to give an entropy contribution of 
the form log \/\R\ = f log L. We also see that this mean- 
field cat state captures the extensive in subsystem size 
number fluctuation while maintaining the ground state 
with definite particle number. 

We can also pin the order parameter to remove the 
anomalous contribution to the entanglement. For the su- 
perfluid, there is now no anomalous entropy although 
there are still extensive number fluctuations. On the 
other hand, in an anti-ferromagnet we may pin the Neel 
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field to point in a particular direction. If we take as a 
mean field state 

|Ncol)=(g)|t)r(g)|4.}r (50) 

reA 

then we see immediately that there is no anomalous en- 
tanglement and the fluctuations of S z are not extensive 
(they are zero). However, the fluctuations of S x and S y 
are still extensive. These simple considerations have been 
considerably developed in Ref. H3- A careful treatment 
including the interactions between the zero mode and the 
sound modes is also expected soon in Ref. [3SL They find 
that the coefficient of the logarithm is somewhat modified 
from the simple minded argument above. 



VI. POTENTIAL VIOLATIONS 

Before concluding, let us point out some potential ways 
that one might violate the scaling forms we have devel- 
oped. A nice example is provided by a certain degener- 
ate spin chain. By fine tuning the strength of the near- 
est neighbor couplings as a function of position along the 
chain, RefAi showed that it was possible to find a ground 
state with entanglement entropy that scaled as the length 
of the interval. This was arranged by adjusting the cou- 
plings so that a real space RG procedure always coupled 
the boundary spin inside the region with a spin outside 
the boundary i.e. no spins formed singlets within the 
interval. This required exponentially decaying couplings 
and is clearly not generic. 

Another route is provided by large ground state de- 
generacy provided we use the completely mixed state. 
Recently such systems have received a lot of attention 
due to results showing non-Fermi liquid behavior in cer- 
tain holographic systems. However, we emphasize that 
nothing is special about the holographic setting, it is but 
one example. A simple condensed matter example with 
"ground state degeneracy" is provided by the spin in- 
coherent Luttinger liquid. This system is a Luttinger 
liquid where the spin energy scale is much less than the 
charge energy scale. At temperatures above the spin en- 
ergy scale, the spin incoherent state emerges where the 
spin degrees of freedom are totally disordered. Such a 
state has an extensive temperature independent entropy, 
but it also cannot be a true ground state. To force the 
state to zero temperature, we must fine tune an infinite 
number of relevant operators (the entire spin Hamilto- 
nian) to zero. We call this an IR incomplete theory since 
it cannot be smoothly connected to zero temperature. 
More generally, we can imagine intermediate scale RG 
fixed points that control the physics over a wide range of 
energies but which cannot be interpreted as ground states 
due to an infinite fine tuning. We know one possibility 
is that such a state may have extensive entropy, but per- 
haps there are other possibilities where the entanglement 
entropy scales like L d ~ 1+S for < S < 1. 



Another setting where violations might occur is in ran- 
dom systems. In one dimension we know that even at in- 
finite randomness fixed points the boundary law for the 
average entanglement entropy is violated is no worse than 
in the conformal case. However, we do not know if infi- 
nite randomness fixed points would violate the boundary 
law in higher dimensions. We expect any finite z random 
fixed point will not, and we suspect that infinite random- 
ness fixed points would not either, but we do not give a 
definite argument at this time. We also note that there 
are considerable subtleties in these systems e.g. typical 
versus average values. Since the entanglement entropy of 
a region has a probability distribution p(S,L), it would 
be interesting to determine if the distribution was a func- 
tion of S/ log (L) only or something more complicated. 
In any event, there are many open issues at such random 
fixed points, the thermodynamics does not obey the sim- 
ple forms we have considered here, and so we do not have 
much else to say about these issues at this time. 

Let us also briefly mention long range interactions. 
If these interactions are due to massless fields with a 
non-singular action within the physical description e.g. 
fluctuating gauge fields or other critical bosonic modes, 
then a proper renormalization group description is pos- 
sible and the entanglement entropy should have no ad- 
ditional anomalous structure. Similarly, so long as the 
long range forces present in the system have such an in- 
terpretation, even if they must be introduced as auxiliary 
fields, we might expect no new anomalies to appear. On 
the other hand, consider the "Id chain" where, in ad- 
dition to nearest neighbor hoppings, every site can hop 
to every other site with the same strength. Calling such 
a system one dimensional is a perversion, but it is an 
extreme form of long range interactions. Clearly such a 
system can be expected to violate the one dimensional 
boundary law more than logarithmically. The task that 
emerges is thus to understand where the crossover point 
is, as a function of the interaction range, to conventional 
one dimensional behavior. One quite interesting situa- 
tion where these considerations are directly applicable 
is momentum space entanglement where the region R 
is some subset of momentum space instead of position 
space. This topic, which has already received some pre- 
liminary attention, deserves its own exposition which we 
will present elsewhere. 

Finally, we note that from the general codimension sec- 
tion our formalism can in some sense encompass states 
that violate the boundary law more seriously than loga- 
rithmically provided that q is effectively less than 1. For 
example, q = roughly describes a state with gapless 
excitations everywhere in some region of finite measure 
in momentum space. This is a lot like the situation with 
ground state degeneracy. Nevertheless, as we have al- 
ready said, we know of no example of a sensible ground 
state with q < 1. 
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VII. DISCUSSION 

We have argued that entanglement entropy and ther- 
mal entropy may be connected via a universal crossover 
function in gapless phases and at critical points. One 
major consequence of this assumption is that local quan- 
tum systems cannot violate the boundary law more than 
logarithmically. However, we hasten to add that should 
our assumptions be violated, we have no objection. In 
particular, possible loopholes escaping our conclusion in- 
clude fine tuning in the Hamiltonian, systems with many 
degenerate ground states, and systems with long range 
interactions. Models showing these characteristics may 
indeed be physically realistic in special cases, neverthe- 
less, we argue that conventional gapless systems, even 
those with critical Fermi surfaces, will not violate the 
boundary law more than logarithmically. Actually calcu- 
lating the entanglement entropy in a model of a critical 
Fermi surface, perhaps in a e expansion, and studying 
in more detail the entanglement properties for q > 1 are 
projects we leave for the future. 

Many conventional quantum critical points that de- 
scribe symmetry breaking transitions, like the 2 + 1 XY 
critical point, fall into the category of conformal field the- 
ories discussed in Sec. II. However, it was recently shown 
in Ref. [3^ that "deconfined" quantum critical points 
have a different entanglement structure due to proximate 
topologically ordered phases. For example, Ref. [H dis- 
cussed the XY* critical point in 2 + 1 dimensions which 
has the same correlation length exponent as the XY tran- 
sition but a different anomalous dimension for the order 
parameter. This arises because the order parameter $ 
is actually composite $ = b 2 and it is the "fractional" 
bosons b that undergo the XY transition. However, an 
important difference is that the b bosons are coupled to 



a Z2 gauge field and hence not all operators in the XY 
theory of the 6s are gauge invariant. Furthermore, there 
must always exist somewhere in the high energy spec- 
trum gapped Z2 vortices. This is all to say that while 
the thermal entropy of XY and XY* are identical, they 
have different crossover functions and different entangle- 
ment properties at zero temperature. 

A potentially profitable generalization of our work here 
would be to include in the scaling formalism the effect 
of relevant operators that move away from the critical 
point. Similarly, it would be interesting to study the 
scaling structure of the full Renyi entropy in more detail. 
Scale invariance fixes the high temperature Renyi depen- 
dence, but the zero temperature entanglement structure 
as a function of Renyi parameter could be rather rich. 
Perhaps our scaling approach could shed some light on 
this structure. The generalization to multiple regions is 
also open and is especially relevant to studies of mutual 
information. 

It would also be quite interesting to develop a vari- 
ational class of density matrices that encode the kind 
of crossover behavior we described here. Since the mu- 
tual information always obeys a boundary law at finite 
temperature^, such states could in principle look like a 
density matrix generalization of tensor network states, 
although presumably the bond dimension would have 
to grow as the temperature was lowered to account for 
systems that violate the boundary law at zero temper- 
ature. Perhaps the new branching MERA approach 40 
could help? One could compute some of these crossover 
functions in field theory to see what sort of universal in- 
formation is easily accessible. We have already done the 
analogous calculations in holographic theories, but there 
it was already known that the entanglement entropy con- 
tains relatively little information due to the large N limit. 
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